Theoretical study of lifetimes and polarizabilities in Ba 
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The 6s~npj (n — 6 — 9) electric-dipole matrix elements and 6s — ndj (n = 5 — 7) electric-quadrupole 
matrix elements in Ba^ are calculated using the relativistic all-order method. The resulting values 
are used to evaluate ground state dipole and quadrupole polarizabilities. In addition, the electric- 
dipole 6pj — 5dji matrix elements and magnetic-dipole 5dc,/2 — 5^3/2 matrix element are calculated 
using the same method in order to determine the lifetimes of the 6pi/2, 6^3/2, 5^3/2, and 6^5/2 
levels. The accuracy of the 6s — 5dj matrix elements is investigated in detail in order to estimate 
the uncertainties in the quadrupole polarizability and 5dj lifetime values. The lifetimes of the 5d 
states in Ba^ are extremely long making precise experiments very difficult. Our final results for 
dipole and quadrupole ground state polarizabilities are asi = 124.15 al and = 4182(34) al, 



respectively. The resulting lifetime values are Tep^ = 7.83 ns, rgp 



P3/2 



;.27 ns, 



81.5(1.2) s. 



and Tsd^^2 ~ 30.3(4) s. The extensive comparison with other theoretical and experimental values is 



carried out for both lifetimes and polarizabilities. 



I. INTRODUCTION 

The atomic properties of Ba"*" ion are of particular 
interest owing to the prospects of studying the parity 
nonconservation (PNC) with a single trapped ion 
Progress on the related spectroscopy with a single Ba+ 
ion is reported in 0, Q , and precision measurements of 
light shifts in a single trapped Ba+ ion have been re- 
ported in [3] . The PNC interactions gives rise to non-zero 
amplitudes for transitions that are otherwise forbidden 
by the parity selection rules, such as 6s— 7s electric-dipole 
transition in Cs. The study of parity nonconservation 
in cesium Q involving both high-precision measure- 
ments and several high-precision calculations provided 
an atomic-physics test of the standard model of the elec- 
troweak interactions and yielded the first measurement 
of the nuclear anapole moment (see 0] for the review 
of study of fundamental symmetries with heavy atoms). 
The analysis of the Cs experiment, which required a cal- 
culation of the nuclear spin-dependent PNC amplitude, 
led to constraints on weak nucleon-nuclcon coupling con- 
stants that are inconsistent with constraints from dee 
inelastic scattering and other nuclear experiments 
More PNC experiments in other atomic systems, such 
as Ba"*", are needed to resolve this issue. The prospects 
for measuringparity violation in Ba+ have been recently 
discussed in Q. 

Ba+ is also of particular interest for developing op- 
tical frequency standard Q and quantum information 
processing [13, [U owing to the extremely long lifetimes 
of 5d states. The accuracy of optical frequency stan- 
dards is limited by the frequency shift in the clock tran- 
sitions caused by the interaction of the ion with exter- 
nal fields. Therefore, knowledge of atomic properties is 
needed for the analysis of the ultimate performance of 
such frequency standard. 

Another motivation for study of Ba+ is an excellent op- 
portunity for tests of theoretical and experimental meth- 
ods, in particular in light of recent measurements of Ba+ 
atomic properties 0, 0, [H, [H, [H, [lB| ■ Ba"*" is a monova- 



lent system allowing for precise theoretical predictions, 
and, in some cases, for evaluation of the theoretical un- 
certainties that do not directly rely on the comparison 
with experiment. It is also an excellent testing case for 
further studies of Ra+ ion, where the correlation correc- 
tions are expected to be larger owing to a larger core. A 
project to measure PNC in a single trapped radium ion 
recently started at the Accelerator Institute (KVI) of the 
University of Groningen (16} . 

In this work, we calculate 6s—npj (n — 6—9), 6pj~5dji 
electric-dipole matrix elements, 6s — ndj (n = 5 — 9) 
electric-quadrupole matrix elements, and 5^5/2 — 5(^3/2 
magnetic-dipole matrix element in Ba^. This set of ma- 
trix elements is needed for accurate calculation of ground 
state dipole and quadrupole polarizabilities and lifetimes 
of the 6pi/2, 6p3/2: 5(i3/2, and 5^5/2 levels. We carefully 
investigate the uncertainty in our values of 6s — 5dj ma- 
trix elements in order to estimate the uncertainties in 
the quadrupole polarizability and the 5dj lifetime values. 
It is particularly important to independently determine 
these uncertainties because of significant inconsistencies 
between different measurements of the 5d^/2 and 5c?5/2 
lifetimes 0, [H [11 [H, [13, [U, HI. There are also 
large discrepancies between experimental determinations 
of the 5d— 6s quadrupole matrix element from the lifetime 
experiments and studies of the Rydberg states of barium 
[12, [H) [l3| • The experimental values of the ground state 
quadrupole polarizability from Refs. [H, [3, differ 
by a factor of two; our value of the quadrupole polar- 
izability is in agreement with Ref. [TJ. We note that 
there are no inconsistencies between the experimental 
lifetimes [l^, [1^, [1^] of the 6pj levels and experimen- 
tal determinations of the electric-dipole ground state po- 
larizability [H [3) [ISl ■ The experimental values of the 
electric-dipole polarizability of the Ba"*" ion in its ground 
state [12, [13, [23 are also in agreement with each other 
and our theoretical value. Our lifetimes of the 6pi/2 and 
6i03/9 levels are in agreement with experimental values 



24|, |25|, l26| within expected accuracy (1%). 
The paper is organized as follows. In Section [Til we 
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give a short description of the method used for the cal- 
culation of the matrix elements. In Section lllli we dis- 
cuss the calculation of the electric-dipole polarizability 
and conduct comparative analysis of the correlation cor- 
rections to the ns — np matrix elements in Ba+, Cs, and 
Ca+. The 6s — Sid quadrupole matrix elements and the 
ground state quadrupole polarizability are discussed in 
Section ITVl and the hfetimes are discussed in Section [V] 
A consistency study of the 5dj lifetime and ground state 
quadrupole polarizability measurements in Ba+ is pre- 
sented in Section Hvl 



II. METHOD 

We calculate the reduced multipole matrix elements 
using the relativistic all-order method [13, [H, 29] which 
is a linearized coupled-cluster method where all single 
and double excitations of the Dirac-Fock wave function 
are included to all orders of perturbation theory. The 
present implementation of the method is suitable for the 
calculation of matrix elements of any one-body operator, 
i.e., the calculations of the El, E2, and Ml matrix el- 
ements are carried out in the same way. We refer the 
reader to the review [2^ and references therein for the 
detailed description of the all-order method. 

Briefly, our starting point is the relativistic no-pair 
Hamiltonian j30| expressed in second quantization as 



ijkl 



(1) 



where a] , aj are single-particle creation and annihilation 
operators, respectively, is the Dirac-Fock (DF) energy 
for the state i, gijki are the two-body Coulomb integrals, 
and : : indicates normal order of the operators with re- 
spect to the closed core. The single-double (SD) all-order 
wave function is written as 
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where is the lowest-order wave function taken to be 
the frozen-core DF wave function of a state v. Indices at 
the beginning of the alphabet, a, b, ■ ■ ■, refer to occupied 
core states, those in the middle of the alphabet m, n, ■ ■ 
refer to excited states, and index v designates the valence 
orbital. The all-order equations for the excitation coeffi- 
cients pma, Pmv, Pmnab, and Pmnva are solvcd itcrativcly 
with a finite basis set, and the correlation energy is used 
as a convergence parameter. The basis set is defined 
in a spherical cavity on non-linear grid and consists of 
single-particle basis states which are linear combinations 
of B-splines [Slj . We use a basis set of 50 splines of order 



9 in a spherical cavity of radius 80 a.u. Such cavity size is 
chosen to accurately represent all orbitals of interest to 
the present study. The resulting excitation coefficients 
Pma, Pmv, Pmnab, and pmnva are uscd to Calculate the 
one-body El, Ml, and E2 matrix elements. 

The SD all-order method yielded results for the pri- 
mary ns — npj El matrix elements of alkali-metal atoms 
that are in agreement with experiment to 0.1%-0.5% 2^. 
We note that while the all-order expression for the matrix 
elements contains 20 terms that are linear or quadratic 
functions of the excitation coefficients, only two terms 
are dominant for all matrix elements considered in this 
work: 



E 

ma 



i^amPwmva + ^maPymwa) 



and 



^ ^ ^ (.^wmPmv H" •^mvPmw) " 



(3) 



(4) 



where Pmnab = Pmnab " Pnmab and z^iv are lowest-order 
matrix elements of the corresponding operator. In the 
case of the electric-quadrupole transitions studied in this 
work, the second term Z^'^^ is overwhelmingly (by an or- 
der of magnitude) larger than any other term. In such 
cases, it was found necessary to include at least partially 
triple excitations into the wave function 
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(5) 

and to correct single excitation coefficient pmv equation 
for the effect of triple excitations [11,133,131, 34]. We have 
conducted such a calculation for the 6s—5dj, 6s—6dj, and 
6s — 7dj electric-quadrupole matrix elements and refer to 
the corresponding results as SDpT values (i.e. including 
all single, double, and partial triple excitations). 

We note that such approach works poorly when terms 
and Z'-"^ are of similar order of magnitude (such as 
all El transition considered here) owing most likely to 
cancellation of high-order corrections to terms and 
Z^'^\ The term Z'"' is not directly corrected for triple 
excitations in the SDpT extension of the method leading 
to consistent treatment of the higher-order correlations 
only when the second term is overwhelmingly dominant. 
We refer the reader to Ref. [35,] for the detailed discussion 
of triple excitations. The results of the matrix element 
calculation are discussed in the following sections. 



III. 



BA+ GROUND STATE DIPOLE 
POLARIZABILITY 



The ground state dipole or quadrupole polarizability 
can be represented as a sum of the valence polarizabil- 
ity ay and the polarizability of the ionic core acore [28| . 
The calculation of the core polarizability assumes allowed 
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excitations to any excited state including the valence 
shell, which requires the introduction of the small counter 
terms a^c to subtract out 1/2 of the contribution corre- 
sponding to the 6s shell excitation \2§\. The core polar- 
izabilities have been calculated in random-phase approx- 
imation (RPA) in Rcf. 36]. The accuracy of the RPA 
values is expected to be on the order of 5% [sij. We cal- 
culated the avc term the in the RPA approximation for 
consistency with acore value. The valence dipole polariz- 
ability for the 6s state of Ba+ is calculated as sum-over- 
states 

1 ( \{6s\\d\\np,/,)\' ^ \{6s\\d\\npy,)\^ \ 
av,Ei = 2^ —f^ ^ + —F^ ^ -(6) 

The sum over the principal quantum number n in Eq. ([5]) 
converges very rapidly and very few first terms have to 
be calculated to high precision. In this work, we use SD 
all-order matrix elements and experimental energies for 
terms with n = 6 — 9 and evaluate the remainder ataii 
in the Dirac-Fock approximation. The contributions to 
the dipole polarizability are summarized in Table IH We 
also list the absolute values of corresponding SD all-order 
reduced electric-dipole matrix elements d. The contribu- 
tion of the terms with n = 6 is overwhelmingly dominant. 
Therefore, the uncertainty in our calculation of the dipole 
polarizability is dominated by the uncertainties in the 
6s — 6pi/2 and 6s — 6^3/2 matrix elements. 

To study the uncertainty in these values, we investigate 
the importance of the contributions from various correla- 
tion correction terms and the overall size of the correla- 
tion correction. The contributions to the 6s — 6pi/2 nia- 
trix element are summarized in Table [III The breakdown 
of the contributions to the 6s — 6^3/2 matrix element is 
essentially the same, and we do not list it here. We also 
give the breakdown of the correlation correction for the 
same transition in Cs and 4s — 'ipi/2 transition in Ca^. 
Cs values are taken from Ref. [33|. Final Ca+ value has 
been published in Ref. ■ As we noted in Section [TTl 
only two terms give large contributions to the correla- 
tion correction. While there are some cancellations in the 
other terms, all them are at least an order of magnitude 
smaller. Unfortunately, there is no straightforward way 
to evaluate the uncertainty in the Z^") term (as we show 
in the later section it can be done for Z'^'^^ ). Therefore, we 
can not make an uncertainty estimate that is independent 
on experimental observations. However, we note that Cs 
6s — 6pj transitions are extremely well studied by a num- 
ber of different experimental approaches (see, for exam- 
ple, [39'| and references therein), and all-order SD data 
are in agreement with Cs experimental values to 0.2%- 
0.4% [231 . The breakdown of terms for Ba+ is slightly 
different than for Cs but is very similar to Ca'*'. As ex- 
pected, the size of the correlations is larger in Ba"*" than 
in Ca"*". Unfortunately, there is only one high-precision 
measurement of the 4pj Ca^ lifetimes [40| that is in sig- 
nificant (2%) disagreement with high-precision theoreti- 
cal results. Similar discrepancies existed for the alkali- 
metal atom measurements done with the same technique 



TABLE I: Contributions to the ground state 6s scalar dipole 
polarizability olbx in Ba^ in units of a%. Comparison with 
experiment and other calculations. The absolute values of 
corresponding SD all-order reduced electric-dipole matrix el- 
ements d (in a.u.) are also given. 



Contribution 


d 




6s - 6pi/2 


3.3357 


40.18 


6s - 6p3/2 


4.7065 


73.82 


6s - 7pi/2 


0.0621 


0.06 


6s - 7p3/2 


0.0868 


0.01 


Oitail 




0.03 






10.61 






-0.51 


Total 




124.15 


Expt. IJA} 




123.88(5) 


Expt. [23] 




125.5(10) 


Theory [41J 




123.07 


Theory [42] 




126.2 


Theory [43] 




124.7 



TABLE II: Contributions of different terms to the Ba^, Ca' 
and Cs ns — npi/2 reduced matrix elements in a.u. 



Contribution 


Ba+ 


Cs [37] 


Ca+ 




6s - 6pi/2 


6s - 6pi/2 


4s - 4pi/2 


DF 


3.891 


5.278 


3.201 




-0.387 


-0.334 


-0.200 


ZM 


-0.209 


-0.485 


-0.120 


Other 


0.041 


0.019 


0.016 


Total 


3.336 


4.478 


2.898 


Correlation 


16.6% 


17.9% 


10.5% 



and later experiments confirmed the theory values. We 
refer the reader to Ref. [38] for more detailed discussion 
of this issue. It would have been very interesting to see 
the 4p lifetimes in Ca"*" remeasured to resolve this issue. 
Based on the similar size of the correlation corrections 
for Cs and Ba"*", we expect similar accuracy of our data 
(on the order of 0.5%). Therefore, the resulting accu- 
racy of our dipole polarizability is expected to be on the 
order of 1%. We find that our value is in excellent agree- 
ment with both experimental values [13, [23j when our 
estimated uncertainty is taken into account. Our results 
are in good agreement with other theoretical calculations 
[41, 42, 43]. We also note that the {6s\d\6p) matrix ele- 
ment has been recently extracted from the K splittings of 
the bound 6snl states in Ref. [ij] , and the resulting value 
{6s\d\6p) = 4.03(12) is in excellent agreement with our 
result (6s]d]6p) = 4.08 (normalized spherical harmonics 
Ci is factored out here for comparison). 
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TABLE III: Absolute values of electric-quadrupole 6s — 5ds/2 and 6s ~ 5d^/2 reduced matrix elements in Ba^ calculated in 
different approximations in a.u. Columns labeled "DF" and "III" are lowest-order Dirac-Fock and third-order MBPT values, 
respectively. The third-order results calculated with maximum number of partial values Imax ~ 6 and Imax = 10 are given to 
illustrate the contribution of the higher partial waves. Breit correction is given separately. The all-order ab initio results are 
given in columns labeled "SD" and "SDpT", respectively; these results include contributions from higher partial waves and 
Breit correction. The corresponding scaled values are listed in columns labeled "SDsc" and "SDpTsc". The calculation of the 
uncertainties of the final values is described in detail in text. 



Transition 


DF 


III Umax = 6) 


III [Irr^ax = 10) 


Breit 


SD 


SDpT 


SD,, 


SDpT.c 


Final 


6s - 5d3/2 


14.76 


11.82 


11.75 


-0.07 


12.42 


12.66 


12.63 


12.59 


12.63(9) 


6s - 5d5/2 


18.38 


14.86 


14.78 


-0.09 


15.55 


15.84 


15.80 


15.76 


15.80(11) 



TABLE IV: Contributions to the ground state 6s quadrupole 
polarizability aE2 in Ba^ and their uncertainties in units of 
uq. The absolute values of corresponding all-order reduced 
electric-quadrupole matrix elements Q (in a.u.) and their 
uncertainties are also given. 



Contribution 


Q 


OtE2 


6s - 5d3/2 


12.63(9) 


1436(20) 


6s - 6d3/2 


16.83(5) 


270(2) 


6s - 7(^3/2 


5.68(5) 


23.7(4) 


6s - 8d3/2 


3.09(6) 


6.3(3) 


6s - 9d3/2 


2.07(4) 


2.7(1) 


6s - 5d5/2 


15.8(1) 


1932(27) 


6s - 6d5/2 


20.30(6) 


392(2) 


6s " 7d5/2 


6.98(6) 


35.7(6) 


6s - 8d5/2 


3.83(8) 


9.6(4) 


6s - 9d5/2 


2.57(5) 


4.192) 


atail 




24(6) 


Otcore 




46(2) 


Total 




4182(34) 


Expt. [14] 




4420(250) 


Expt. [12] 




2462(361) 


Expt. [23] 




2050(100) 



IV. BA+ GROUND STATE QUADRUPOLE 
POLARIZABILITY 

The valence part of the quadrupole polarizability is 
given in the sum-over-states approach by 

1 ^ / \{Qs\\Q\\nd^/2)? ^ \{&s\\Q\\nd,/^)\^ \ 

The RPA core value [s^ is 46aQ, and the a^c term is neg- 
ligible. The terms containing the 6s — 5c?3/2 and — 
matrix elements give overwhelmingly dominant contri- 
bution to the total values. Therefore, we study these 
transitions in more detail and evaluate their uncertain- 
ties. Unlike the case of the El transitions considered 
earlier, Z^'^'^ term contributes over 90% of the total cor- 
relation correction. Therefore, we carried out the cal- 
culation using both SD and SDpT approaches described 



in Section [TTl We also carried out semi-empirical scal- 
ing in both approximations by multiplying single exci- 
tation coefficients Pmv by the ratio of the "experimen- 
tal" and corresponding (SD or SDpT) correlation ener- 
gies [13. The "experimental" correlation energies are 
determined as the difference of the total experimental 
energy and the DF lowest-order values. The calculation 
of the matrix elements is then repeated with the modi- 
fied excitation coefficients. The accuracy of such scaling 
procedure for the similar cases was discussed in detail 
in Refs. [H, [s^, [13] ■ The reasoning for such a scaling 
procedure in third-order perturbation theory (scaling of 
the self-energy operator) has been discussed in Ref. [45| . 
We list SD, SDpT, and the corresponding scaled results 
(labeled "SD^^" and "SDpT^^") in Tablc[nD The lowest- 
order DF results are listed to illustrate the size of the cor- 
relation corrections. We demonstrate the size of the two 
other corrections, contribution of the higher partial waves 
and Breit correction, in the same table. The first correc- 
tion results from the truncation of the partial waves in all 
sums in all-order calculation at Imax — 6. All-order calcu- 
lation with higher number of partial waves is unpractical. 
Therefore, we carry out the third-order MBPT calcula- 
tion (following Ref. ^]) including all partial waves up 
to Imax = 6 and Imax = 10 and take the difference of 
these two values to be the contribution of the omitted 
partial waves that we add to ah initio all-order results. 
We verified that the contribution of the Z = 9 — 10 partial 
waves is very small justifying the omission of contribu- 
tions from I > 10. The Breit correction is calculated as 
the difference of the third-order results with two different 
basis sets. The second basis set is generated with taking 
into account one-body part of the Breit interaction. We 
note that scaled values should not be corrected for ei- 
ther partial wave truncation error or Breit interaction to 
avoid possible double-counting of the same effects. We 
take SDsc values as our final results. The uncertainty 
of the final values is calculated as follows: the uncer- 
tainty in the Z^'^' term is evaluated as the spread of the 
most high-precision values (SDsc, ab initio SDpT, and 
SDpTsc), the remaining theoretical uncertainty in the 
Coulomb correlation correction is taken to be the same 
as the uncertainty in the dominant Z^"^^ term. We assume 
100% uncertainties in the contributions of the higher par- 
tial waves and Breit correction. The final uncertainty of 
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the 6s — 5dj matrix elements (0.7%) is obtained by adding 
these four uncertainties in quadrature. We note that this 
procedure for the uncertainty evaluation does not rely on 
the experimental values with the exception of the exper- 
imental energies used for scaling. 

The contributions to the ground state quadrupole po- 
larizability are given in Table [TVl While the n = 5 term 
is dominant, the contributions of the few next terms are 
substantial. Therefore, we carry out SD, SDpT, and 
both scaled calculations for the 6s — 6dj and 6s — 7dj 
matrix elements as well and repeat the uncertainty anal- 
ysis described above (we omit Breit and higher-partial 
wave corrections here since such precise evaluation of the 
uncertainties is not needed for these transitions). The 



6s 



8dj and 6s 



9c?,- matrix elements are calculated in 



third-order MBPT, and their accuracy is taken to be 2% 
based on the comparison of the third-order and all-order 
values of the 6s — 7dj matrix elements. The remainder is 
evaluated in the DF approximation and reduced by 23% 
based on the comparison of the DF and third-order data 
for 6s — 8dj and 6s — 9dj matrix elements. Its accuracy 
is correspondingly taken to be 23%. 

Our recommended value for the ground state 
quadrupole polarizability is in agreement within the cor- 
responding uncertainties with the most recent experi- 
mental work p^ . However, our value for the contribu- 
tion of the 6s — 5dj transitions to the quadrupole po- 
larizability [3368(34)] differs by about a factor of 2 from 
the experimental values [l^, [3, [l!] obtained based on 
the nonadiabatic effects on the Rydberg fine-structure 
intervals. This issue and the discrepancies in the ex- 
perimental values of the quadrupole polarizabilities are 
addressed in detail in Ref. [14]. We note that these ex- 
perimental values of the 6s — 5dj contributions to the 
quadrupole polarizabilities (1524(8) 0] and 1562(93) 
[l3| in the two most recent studies) are significantly in- 
consistent with all high-precision calculations of the 5c?,- 



lifetimes [15|, IM |47|, |48|, |49|, |50| carried out by differ- 
ent methods as well as with all experimental lifetime 
measurements (also carried out by different techniques) 
[H, M, M, Hi, H iH, m. For comparison, the value 
1562(93) obtained from the (6s|r2|5c?) = 9.76(29) matrix 
element that was extracted from the K splittings of the 
bound 6sn^ states in Ref. [l^ corresponds to the lifetime 
T5d3/2 = 170(10) s that is a factor of 2 longer than all 
other values. We discuss the lifetimes of the 5^3/2 and 
5c?5/2 levels in the next section. 



V. LIFETIMES 

The lifetime of a state a is calculated as Tq = 
iJ2b<a ^ab)~^- The El, E2, and Ml transition rates Aab 
are given by [5l| : 



2.02613 X 10^^ S 



El 



A3 



2.7 a + 1 



(8) 



TABLE V: Lifetimes of the 



and 5dj states in Ba 



parison with experiment and other theory, 
the 6pj states are given in ns, 
states are given in s. 



com- 

The lifetimes of 
and the lifetimes of the 5d, 





^6pi/2 (ns) 


r6p3/2 (ns) 




^5d5/2 (S) 


Present 


7.83 


6.27 


81.5(1.2) 


30.3(4) 


Theory [46. 471 


7.99 


6.39 


83.7 


30.8 


Theory [48] 


7.89 


6.30 


81.5 


30.3 


Theory [49] 


7.92 


6.31 


81.4 


36.5 


Theory [50| 






80.1(7) 


29.9(3) 


Theory [151 






82.0 


31.6 


Expt. [241 




6.312(16) 






Expt. [25] 


7.92(8) 








Expt. [26] 


7.90(10) 


6.32(10) 






Expt. [17] 






17.5(4) 




Expt. [18] 






48(6) 




Expt. [19| 






79.8(4.6) 




Expt. \15} 






89.4(15.6) 


32.0(4.6) 


Expt. [20| 








47.0(16) 


Expt. [21] 








32.0(5) 


Expt. [22] 








34.5(3.5) 



aE2 

^ab 



1.11995 X 10 



18 



Se2 



A 



Ml 
b 



A5 2ja + 1 

2.69735 X 1013 Sm 



A3 



2ia + 1 



(9) 
(10) 



respectively, where A is the wavelength of the transition 
in A and S is the line strength. In this work, we calcu- 
lated the lifetimes of the Qpi/2, 6^3/2, 5c?3/2, and 5^5/2 
levels in Ba"*" . The results are compared with experimen- 
tal and other theoretical values in Table |Vl Since the 6p 
levels are above 5c? levels in Ba"*", we also needed to cal- 
culate the SD all-order reduced matrix elements for the 
6p — 5c? El transitions, and our results (in atomic units) 
are d(6pi/2 - 5c?3/2) = 3.034, c?(6p3/2 - 5c?3/2) = 1.325, 
and c?(6p3/2 — 5c?5/2) = 4.080. These values include con- 
tributions from the higher partial waves (0.6%) and 0.1%- 
0.2% Breit correction. The correlation corrections to 
these transitions are similar to the ones for the 6s — Qpj 
transitions. Therefore, similar (on the order of 0.5%) 
accuracy is expected for these matrix elements. The 
6s — &pj transitions contribute about 73% to the respec- 
tive X]b<a -^a-b totals for the &pj lifetimes. Based on our 
evaluation of the uncertainty in these matrix elements 
discussed in Section IIIIl we expect present 6p lifetime 
values to be accurate to about 1%. Our results are in 
excellent agreement with other recent theoretical [3 El] 
and experimental [U, [2^ [2^ values. The calculation 
of Refs. [H, 113] is a third-order MBPT calculation that 
omits higher-order corrections included in the present 
calculation, slightly different values are expected. 

Only one transition contributes to the 5c?3/2 lifetime: 
6s — 5p3/2 E2 transition (the contribution of the 6s — 5c?3/2 
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Ml transition is negligible). In the case of the 5^5/2 life- 
time, Ml 5(^5/2 — 5( ^3/9 transition has to be included as 
pointed out in [13, l48l . Iso! ] . Our SD all-order value for 
this transition (in a.u.) is 1.5493. The correlation cor- 
rection contribution is very small, and the lowest order 
gives essentially the same value, 1.5489. The Ml transi- 
tion contributes 17% to the X]6<a total for the 5^5/2 
level. 

We compare our final results for the 5d^/7 and 5cfc^ 
lifetimes with experimental R|, [ll Jlljm 0, HH [13 
and other theoretical [IlliEllaJlaltils^values in Ta- 
ble|Vl We note that calculation [49| omitted 5d^/2 — 5(1^/2 
Ml contribution to the 5c?5/2 lifetime leading to a higher 
value, as noted in later work 150]. Our results are in 
agreement with other theoretical calculations, most re- 
cent values from measured in a beam-laser experi- 
ment performed at the ion storage ring CRYRING, as 
well as experimental values from [l^, l2ll [2^ . 



elements; and 5(^5/2 — M3/2 magnetic-dipole matrix el- 
ement. These values are used to evaluate lifetimes of the 
6pi/2j 6p3/2, 5(^3/2, and 5^5/2 levels as well as dipole and 
quadrupole ground state polarizabilities. Extensive com- 
parison with other theoretical and experimental values 
is carried out. The present values of the dipole polariz- 
ability and 6pj lifetimes are in excellent agreement with 
experimental values. We estimated the uncertainty of our 
theoretical values for these properties to be on the order 
of 1%. Our recommended value of the quadrupole ground 
state polarizability aE2 — 4182(34)aQ is in agreement 
with the most recent experimental work [14']. Our rec- 
ommended values for the 5dj hfetimes rsdg/a — 81.5(1.2) s 
and T5d^/2 ^ 30.3(4) s are in agreement with other the- 
oretical calculations, most recent values from [l5| mea- 
sured in a beam-laser experiment performed at the ion 
stora ge ring CRYRING, as well as experimental values 
from]ii,[2l[23]. 



VI. CONCLUSION 

In conclusion, we carried out the relativistic all-order 
calculations of Ba"*" 6s — npj (n = 6 — 9), 6pi/2 — 5^3/2, 
6p3/2 —5^5/2, and 6^3/2 — 5d5/2 electric-dipole matrix el- 
ements; 6s — 5^3/2, 6s — 5(i5/2, 6s — 6^3/2, 6s — 6^5/2, 
6s — 7c?3/2, and 6s — 7^5/2 electric-quadrupole matrix 
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